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Abstract 

We show that quantum probabihties can be derived from statistical 
mechanics of classical fields. We consider Brownian motion in the space 
of fields and show that such a random field interacting with threshold 
type detectors produces clicks at random moments of time. And the cor- 
responding probability distribution can be approximately described by 
the formalism of quantum mechanics. Hence, probabilities in quantum 
mechanics and classical statistical mechanics differ not so much as it is 
typically claimed. The temporal structure of the "prequantum random 
field" (which is the L2-valued Wiener process) plays the crucial role. Mo- 
ments of detector's clicks are mathematically described as hitting times 
which are actively used in classical theory of stochastic processes. Born's 
rule appears as an approximate rule. In principle, the difference between 
the "precise detection probability rule" derived in this paper and the con- 
ventional Born's rule can be tested experimentally. In our model the 
presence of the random gain in detectors playes a crucial role. We also 
stress the role of the detection threshold. It is not merely a technicality, 
but the fundamental element of the model. 

keywords: classical statistical mechanics, infinite-dimensional state space, 
random fields, quantum probability of detection, derivation of Born's rule, 
threshold detectors, asymptotics of error function, distribution of hitting times, 
Brownian motion in the space of fields 

1 Introduction 

In this paper we present a novel application of statistical mechanics with the 
infinite-dimensional state space, the space of classical fields. We study the 
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problem of interaction of a random field, e.g., the electromagnetic field, with 
a detector of the thershold type. And we found that, for the natural random 
field, the Brownian motion in the space of fields (the L2-valued Wiener process), 
statistics of clicks produced by such a detector can be approximately described 
by the Born's rule. In quantum mechanics (QM) the latter was postulated [T]. 
Whether it is possible to derive this rule from some natural physical principles 
is still the subject of intensive debates^ Thus we show that classical statistical 
mechanics of fields (theory of random fields [B]-[TD]) provides a possibility to de- 
rive the fundamental law of QM connecting theory with experimental statistics. 

We state again that the basic quantum rule is derived as an approximate 
rule. (Opposite to the conventional quantum mechanical viewpoint; by the latter 
this rule is a "precise rule".jl The small parameter of the model (determining 
the magnitude of deviation of the Born's rule from the "precise rule" for the 
threshold type detectors interacting with the L2-valued Wiener process) is the 
quantity 

e^%^«l, (1) 

where £d is the detection threshold (having the physical dimension of energy) 
and Impulse is the average energy of pulses emitted (randomly) by a source. Thus 
the theory is about random fields of low energy interacting with detectors with 
sufficiently high threshold. We derived the "precise formula" for the detec- 
tion probability for the threshold type detectors. The basic idea behind this 
derivation is embedding of the problem of detection moments for the threshold 
detection into well developed theory of hitting times |12j. (The latter is widely 
used in e.g. statistical radio-physics [13].) This theory gives us the probability 
distribution of detection times. The final (precise) formula, see ((38)) . is quite 
complicated analytically. At the same time it is very simple from the view- 
point of numerical computation (summation of a series with quickly decreasing 
terms). In principle, this formula can be tested experimentally^ 

This paper can be considered as development of prequantum classical statisti- 
cal field theory (PCSFT) [ll]-[23]. PCSFT provided representation of quantum 
averages and probabilities as averages and correlations with repspect to en- 
sembles of classical fields; including entangled systems (cf. IIofer|4],de la Pena 
and Ceto [24], Casado et al. 25 , Boyer [26], Cole [22], Roychoudhuri [28], 
Adenierf29J). So, QM and classical statistical mechanics became essentially 
closer than it was commonly believed. All quantum probabilistic quantities are 



^ "The conclusion seems to be that no generally accepted derivation of the Born rule has 
been given to date, but this does not imply that such a derivation is impossible in principle." 

[2],cf. [3], g], m- 

^We remark that recently the quantum optics group of G. Weihs performed a series of 
experiments testing the validity of the Born's rule in the triple-slit experiment . Whether 
the observed deviation from the Born's rule is of the fundamental nature or it is a consequence 
of experimental technicalities such as e.g. nonlinearity of detectors is still an open question. 

•^The precise formula for the threshold detection probability contains the probability dis- 
tribution of the detector's gain. Determination of this distribution is a complicated problem 
depending essentially on types of threshold detectors [?]. We shall discuss this problem in 
section IT2I 
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represented in the classical probabilistic framework, cf. Manko et al. |53|-|38]. 
Hess et al. [39], gl], De Raedt et al. [40], Garcia et al. [42]. PCSFT reproduced 
these quantities as averages of intensities of continuous random fields. However, 
quantum experimental statistics is based on discrete events, clicks of detectors. 
To come closer to experiment, PCSFT should be completed by measurement 
theory describing the transition from continuous random intensities to discrete 
random events of detection. 

The first model of the thershold detection (TSD) of continuous random fields 
serving PCSFT was proposed in [33]. This model is mathematically compli- 
cated. The "prequantum random fields" reproducing quantum statistics are 
very singlular (of the white noise type). The state space of this model is the 
space of Schwartz tempered distributions. In [44], it was observed that the 
same result can be achieved by using essentially simpler model, namely, Brown- 
ian motion (Wiener process). However, in j44) from the very beginning a rough 
asymptotics of the detection probabilities was used; hence, directly the Horn's 
rule (an approxiamte formula for detection probabilities) was derived. In the 
present paper we use deeper results from theory of classical random fields (for 
hitting times) and obtain the "precise rule" for threshold detection which can 
be compared with the Horn's rule, the approximate rule. Moreover, in |44) only 
internal degrees of freedom such as e.g. polarization were considered. In the 
present paper by starting with internal degrees of freedom we proceed to the 
really physical case of random fields on physical space. (The latter model is 
classical statistical mechanics with infinite-dimensional state space.) 

As was already stressed, our approach is based on embedding of the de- 
tection problem into theory of random hitting times. As a consequence, we 
shall use the mathematical apparatus of classical probability and mathematical 
statistics. The error function of the Gaussian distribution and compllementary 
error function (its asymptotic expansion) are our main mathematical tools, see, 
e.g., [45] for introduction and f46; for applications. 

2 Threshold detectors 

The detection procedure used in optical experiments is based on two steps 
threshold passing. The first threshold is so to say the "hardware threshold". 
This is a part of the solid state physical structure of the detector, we call it 
the fundamental threshold. The second threshold is in some sense the "software 
threshold" - the discrimination threshold. This threshold is set by experiments, 
in the case of photo-multipliers and tungsten-based superconducting transition- 
edge sensors (W-TESs), or directly in the process of fabrication, in the case 
of silicon-avalanche-photodiodes [47]. In this paper we are basically interested 
in the discrimination threshold. However, we start with a brief discussion on 
fundamental thresholds. Hence, everywhere besides section \27l\ the "detection 
threshold" is just the discrimination threshold. 
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2.1 Fundamental thresholds of detectors 

There is a difference in tliresliold's structure between photo-multipliers, as tfie 
one used in [SD], [ST], [52], and more recent experiments (e.g., [S3]-[57]) based 
on silicon- avalanche-photodiodes. However, detectors of botfi types are based on 
tfiresfiold processes. Ttie total tfiresfiold is combined of a few different tliresh- 
olds corresponding to different stages of tlie detection process. First, we de- 
scribe fundamental thresholds corresponding to plioto-multipliers and silicon- 
avalanche-pfiotodiodes. In photo-multipliers the fundamental threshold is given 
by the work function, i.e., ffund coincides with the work function. In silicon- 
avalanche-photodiodes the fundamental threshold is given by the band gap, i.e., 
iffund coincides with the band gap. 

2.2 Discrimination threshold 

A lot of noise is involved in the process of detection, e.g. [47]-[49], An important 
source of noise is the multiplication process. For example, in photo-multipliers, 
once an electron has been extracted from the metal, it is accelerated in vacuum 
by an electric field until its kinetic energy is enough to extract other bounds 
electrons (secondary emission) when striking the surface of another metal surface 
(dynode), which will in turn be used accelerated onto other dynodes to free 
more and more electrons, until that flow of electrons becomes measurable as 
anodic current. The form and energy of output spikes vary significantly from 
one liberated energy carrier to another. Typically it is assumed that the gain is 
given by 

G = (2) 

where ^ is the multiplication factor of a single dynode, a is the fraction of 
photoelectrons collected by the multiplier structure, and N is the number of 
stages in the photomultiplier. In the most simple model, ^ can be assumed to 
follow a Poisson distribution about the average yield for each dynode, so that 
the gain is a compound Poisson process over N identical stages. However, ex- 
perimental measurements of the single photoelectron pulse height spectra from 
photo-multipliers exhibit a distribution with larger relative variance than pre- 
dicted by the Poisson model (and in some case with a decreasing exponential 
distribution instead of a peaked one) , and there is thus no universal description 
of multiplication statistics [47] . 

Besides of noise produced by detectors, noise responsible for so called dark 
counts plays an important role. This noise of the random background is in- 
escapable. In average spikes corresponding to signal detections differ in the 
amplitude from noise generated spikes. This fact provides a possibility to filter 
noise generated spikes by using a discrimination threshold, denote the later by 
the symbol £d- The selection of this threshold is a delicate procedure. By se- 
lecting too low £d experimenter would count to many noise generated spikes, in 
particular, dark counts. By selecting it too high experimenter would discard too 
many spikes generated by the signal. In both cases quantum statistics would be 
essentially disturbed. 
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During a last few years, a few leading groups in quantum optics started 
to use actively with Tungsten-based Superconducting Transition-Edge Sensors 
(W-TESs) - the ultra-sensitive microcalorimeters (in particular, in attempts to 
perform the EPR-Bell experiment with highly efficient detectors). Suprisingly 
usage of such modern detectors is also impossible without seting a proper dis- 
crimination threshold, see, e.g., [15], [3^. for reviews. I guess that (since noise 
is everywhere) it is impossible to elaborate a photon detection procedure which 
is not of the threshold type. 

3 Threshold detection scheme 

3.1 Detection moment as hitting time 

We consider a threshold type detector with the threshold £d- It interacts with 
a random field 0(s,aj), where s is time and is a chance parameter describing 
randomness. For a moment, we consider the C-valued random field (complex 
stochastic process) . In section 14.11 we shall consider random fields valued in 
finite-dimensional complex space H. This correspond to detection of internal 
degrees of freedom such as e.g. polarization. We stress that the real physical 
situation corresponds to random fields with infinite-dimensional state space, 
e.g., H = L2(R'^), the space of complex valued fields : — C (or — )• C*"'), 
see section ?? 

The energy of the field is given by £{s,uj) — |0(s,w)p (hence, the random 
field has the physical dimension ~ ^energy) . (In section 5 we shall also consider 
spatial variables. In this case the field has the physical dimension of the energy 
density, i.e., energy per volume.) A threshold detector clicks at the first moment 
of time t{uj) when signal's energy £ multiplied by the gain g, see section [2.21 
exeeds the threshold: 

g£iTiLu),Lu)>£d. (3) 

In the mathematical model the detection moment is defined as the first hitting 
time; 

t(w) = inf{.s > : g£{s, uj) > £d}. (4) 

3.2 PDF of the detection moment 

We proceed under the following basic assumption. After arriving to a threshold 
type detector a classical signal (random field) behaves inside this detector as 
the (complex) Brownian motion, i.e., the (j){s,uj) is simply the Wiener process, 
the Gaussian process having zero average at any moment of time 

E<j){s,uj) =Q. (5) 

and the covariance function 

E4){si,uj)cl){s2,uj) = min(si,S2)cr^; (6) 
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in particular, we can find average of its energy 



E£{s,Lu) =a^s. (7) 

From this equation, we see that the coefficient has the physical dimension of 
power. We are interested in the probability distribution of the moments of the 
fd-threshold detection for the energy of the Brownian motion. Since moments 
of detection are defined formally as hitting times, we can apply theory of hitting 
times [12], [8]. Consider 

Ta{co) = mi{s > : £{s, uj) > a^} = inf{s > : \(f>{s, uj)\ > a}. (8) 

Its probability distribution function (PDF) is given by the complicated expres- 
sion, see, e.g., Shyryaev [5]: 

P(r. < At) = [l - *(^^^^)] , (9) 

fc=o va I\t 

where $(a;) = e"^^ /"^du is the PDF of the standard Gaussian distri- 

bution. To show close relation to classical statistics, we also consider the error 
function [45] 

erf(x) - ^ / e~^'du 
V"" Jo 

and the complementary error function [l^ 

2 f°° _ 2 
erfc(x) = 1 — erf(x) = —= / e " du. 

Now we select a ~ ^ — and set r = Tq. We obtain: 



Piricu) < At) ^ 2^(-l)'=crfc((l + 2k) J^-^). (10) 



fc=0 



3.3 Asymptotics for the detection moment 

In coming considerations. At is the average duration of the interaction of a sig- 
nal with a threshold detector. The experimental scheme can be described in the 
following way. There is a source of random pulses of e.g. classical electromag- 
netic field. (Such pulses can be identified with wave packets used in the quantum 
formalism.) Each pulse propagates in space (we shall discuss its dynamics) and 
finally arrives to a detector. The aforementioned temporal parameter At is the 
(average) duration of interaction of an input pulse with a detector. The quantity 

f pulse = cr^Ai (11) 
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is the average energy of emitted pulses0 We shall proceed under the basic 
assumption that this energy is essentially less than the threshold, i.e., 

e=%^«l. (12) 

This is a realistic assumption, since the discrimination threshold is set for the 
amplified signals from the detector and the gain producing this amplification is 
very large. Under this assumption the probability distribution P{t{uj) < At) 
can be approximated by the first term in the series ([1]): 

P(T(tj) < At) w 2erfc(l/V2eg). (13) 

In further consideration we shall use the asymptotic behavior of the comple- 
mentary error function [IS] for a; — t- oo : 



erfc(x) 



(2x2) 

n=l ^ ' 



i+f(-i)n ^-3-5--(^"-^) .£r::f(-i)n(^ 



(14) 

To proceed towards derivation of the Born's rule, we are fine with the rough 
approximation 

erfc(x) « (15) 

X^TT 

(So, in our approach the Born's rule provides a rather rough approximation of 
detection probabilities.) Hence, 



3.4 Asymptotics in the presence of a random gain 

As was pointed in section 12. 2[ the gain is by itself a random variable [] . Thus 

t{uj) =: inf{.s > : g{u)£{s,u) > Sa}- (17) 
Let us introduce the random variable rj by setting 

9 = I/'?'- 

Thus in the presence of a random gain the detection moment is given by t{uj)) = 
inf{s > : £{s,uj) > ri'^{u!)£d}- Denote the probability distribution of by P,,. 
We suppose that it has the density, denote the latter by the symbol Pr)(A). To 
find the probability distribution of t, we use the formula of total probability: 



P{t <At) = j d\pr,{\)P{T < Atlij = A), 



(18) 



fact, this is the average energy which is transmitted to a detector by a pulse in the 
process of interaction. We proceed under the assumption that there are no losses and a 
detector "eats" (in average) the total energy of the emitted pulse. 
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where the conditional probability has been already found: 



P(.<A.|,^A,.2^55|eM-^). (19) 

Finally, we have 

£d J |A| ^2^£i^ 
3.5 The average number of clicks 

Wc now recall the meaning of the time intcirval At : the average duration of the 
interaction of a signal with a detector. Consider now a long run of experimental 
trials, emissions and detections, of the duration T. Then the average number of 
clicks can be found as 



where /,,(A) = We now remark that 



exp{-^} 



V2tt€ 

is the (^-sequence for e ^ : lime_>.o -De(A) = (5(A). Hence, 

4 Detection probability 

4.1 Derivation of the quantum (Born's) rule for the de- 
tection probability 

Consider now a random signal (/)(s, lo) valued in the m- dimensional complex 
Hllbert space H, where m can be equal to infinity. Let {ej) be an orthonormal 
basis in H. The vector- valued (classical) signal 0(s,a;) can be expanded with 
respect to this basis 

(f>{s,uj) = ^(pj{s,uj)ej, (f>j{s,uj) = {ej,(f>{s,uj)). (23) 
j 

This mathematical operation is physically realized as spliting of the signal 
(f>{s,Lu) into components (j)j{s,uj). These components can be processed through 
mutually disjoint channels, j = 1, 2, ...m. We now assume that there is a thresh- 
old detector in each channel, Di, Dm- We also assume that all detectors have 
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the same threshold f ^ > and the same probabihty distribution of the gain, 
with the density Pg{X)- 

Suppose now that <j}{s,ijj) is the valued Brownian motion (the Wiener 
process in H) . This process is determined by the covariance operator B : H ^ 
H. Any covariance operator is Hermitian, positive, and the trace-class and vice 
versa. The complex Wiener process is characterized by the Hermitian covariance 
operator. (We remark that complex-valued random signals are widely used in 
e.g. radio-physics.) We have, for y & H, 

E{y,cp{s,Lo)) 

and, for yj e H,j — 1, 2, 

E{yi,(t){si,uj)){(t){s2,uj),y2) = mm{si, S2){Byi,y2). 

The latter is the covariance function of the stochastic process; in the operator 
form: B{si,S2) — mm{si, S2)B. We note that the dispersion of the iJ-valued 
Wiener process (at the instant of time s) is given by 

= £;||0(s,a;)|p sTrB. (24) 

This is the average energy of this random signal at the instant of time s. Hence, 
the quantity 

= ^ = TrB (25) 

s 

has the physical dimension of power; this is the average power of the signal of 
the Brownian motion type (it does not depend on time). 

We also remark that by normalization of the covariance function for the fixed 
s by the dispersion we obtain the operator, 

p = S/TrB, (26) 

which formally has all properties of the density operator used in quantum theory 
to represent quantum states. Its matrix elements have the form 

p,, = 6,j/e2. (27) 

These are dimensionless quantities. The relation plays a fundamental role 
in our approach |14j-|23|: each classical random process generates a quantum 
state (in general mixed) which is given by the normalized covariance operator 
of the process. One can proceed the other way around as well: each density 
operator determines a class of classical random processes [58] . 

We now consider components of the random signal 0(s,cj) and their corre- 
lations: 



E(j)i{si,uj)(l)i{s2,uj) ^ mm{si,S2){Bei,ej) ^ min{si, S2)hj . (28) 
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In particular, 

cr|(s) = E£j{s,uj) = E\(l)j{s,uj)f s bjj. (29) 

This is the average energy of the jth component at the instant of time s. We 
also consider its average power: 

a] = b,,. (30) 

We remark that the average power of the total random signal is equal to the 
sum of powers of its components: 



^-E-l (31) 



Consider now a run of experiment of the duration T. The average number 
of clicks for the jth detector can be approximately expressed as 

5!^. (3.) 

This asymptotics is valid under the assumption (jlip. Here 

is the average energy of emitted pulses. (We state again that we consider the 
process of detection in the absence of losses, cf. footnote 3.) 
The total number of clicks is (again approximately) given by 

N^Y.^,^'^!^. (33) 

3 

Hence, for the detector Dj, the probability of detection can be expressed as 

This is, in fact, the Bom's rule for the quantum state p and the projection oper- 
ator Cj = |ej)(ej| on the vector ej : Hence, for the detector Dj, the probability 
of detection can be expressed as 

P,=TrpCj. (35) 
4.2 Ercf-rule for detection probability 



We state again that the expression ()22|) was based on an approximation of the 
PDF of the detection moment for the detection threshold £d» i^^puisc- 

By using the precise formula (jlOp and by taking into account the presence 
of the random gain we obtain 
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^£(-1)^/ ciAp,(A)erfc((l + 2k) /^). (36) 



k=0 V J 

The total number of clicks in all detectors is the sum of Nj 



N-^^t^i-^)"J dAp,(A)^erfc((l + 2k) /^). (37) 

k=0 j V j 

Hence, the probability of a click in the jth detector is given by sufficiently 
complex formula (generalized Born's rule): 



EZoi-^r I dAp,(A)erfc((l + 2k)y^ 



„ _„?At, 

EZoi-^)' I d\Pn{\) erfc((l + 2k)y^' 



We can easily rewrite this formula by using the original probablity distri- 
bution of the gain g = g{uj) and not of the random variable ri{uj). (We state 
again that the latter random variable is related to the gain by the relation 
g{ui) = l/rf'{ijj).) We remark that densities of the probability distributions for 
g and "q (for nonnegative 77) are connected by the rule: 

1^ j d\p,{\)P{T, < M\g ^ \) 



We have 



gf:(-l)^/ .Ap,(A)erfc((l + 2k) /^). (40) 



fe=o V J 

The total number of clicks in all detectors is the sum of Nj 



2At 



(41) 



Hence, the probability of a click in the jth detector is given by sufficiently 
complex formula (generalized Born's rule): 



Er=o(-l)'/rfV.(A)erfc((l + 2k)y=lf^^ 
' Er=o(-l)'=/rfAp,(A)E, erfc((l + 2k)y=|^^ ' 



As we have seen, it is not easy to derive the canonical Born's rule from 
(|38p . Although this expression is very complex from the analytic viewpoint, 
it can be easy and quickly calculated by computer. It would be interesting to 
compare the result which can be obtained on the basis of (1381) with experimental 
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data. TSD predicts that the experimental statistical data should match with 
this generalaized Bern's rule better than with the canonical Bern's rule. 

We also remark that, for sufficiently large values of f^, the series inside 
expressions (pS)) , coverge very quickly. Therefore by using only the first 
terms of these series we obtain very good approximation. Hence the following 
approximate Ercf-expressions for the detection probabilities can be used 



/dA/9^(A)2.crfc 



2a'^At , 



or directly in the gain variable 



Pj « \ / N ■ (44) 

/dAp,(A)E,crfc(./=f^ 



We remark that, although the gain variable has the straightforward physical 
meaning, and hence the representation (j44l) is better from the physical view- 
point, the 77-representation is better from the mathematical viewpoint: the 
integrals in (j43]) are Gaussian. The same can be said about complete Ercf- 
representations of the detection probabilities, cf. (j42|) with ((38|). 
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